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Abstract 

Cosmological perturbations generated quantum-mechanically (as a particu- 
lar case, during inflation) possess statistical properties of squeezed quantum 
states. The power spectra of the perturbations are modulated and the angular 
distribution of the produced temperature fluctuations of the CMBR is quite 
specific. An exact formula is derived for the angular correlation function of 
the temperature fluctuations caused by squeezed gravitational waves. The 
predicted angular pattern can, in principle, be revealed by the COBE-type 
observations. 
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The recent discovery by COBE [1] of the angular variations of CMBR makes it neces- 
sary to analyze in greater detail the observational consequences of the quantum-mechanical 
generation of cosmological perturbations. The underlying physical reason for the generating 
process is the parametric (superadiabatic) amplification of classical perturbations and the 
associated quantum- mechanical particle pair creation in the variable gravitational field of the 
homogeneous isotropic Universe. As a result of the parametric coupling between the quan- 
tized perturbations and the variable classical "pump" field, the initial vacuum state of the 
perturbations evolves (in the Schrodinger picture) into a strongly squeezed vacuum state 
possessing very specific statistical properties. The generated fluctuations can be viewed, 
classically, as a stochastic collection of standing waves. The mechanism itself and its main 
results concerning squeezing are valid for gravitational waves and progenitors of density per- 
turbations [2,3]. A particular variable gravitational field, that may be responsible for the 
amplification process, is provided by one or another type of the inflationary expansion. It 
is often stated that inflation generates "Gaussian perturbations with randomly distributed 
phases". However, this is not the case: the phases of all modes of perturbations are es- 
sentially constant and fixed [3] which leads to standing waves, modulated spectra of the 
generated perturbations and a specific angular distribution of the temperature fluctuations 
of CMBR over the sky, as will be shown below. 

In this paper we will analyze, mostly, gravitational waves. For our purposes it is sufficient 
to consider perturbations in a spatially-flat FLRW universe 
ds^ = a^{rj){drf — dx^ — dy^ — dz^) where airj) is the cosmological scale factor. 

The quantum-mechanical operator for the gravitational- wave field can be written in the 
general form 

2 

/oo ^ 
rf'nEp|,(n)K(,7)e^"- + a^^{r^)e-'-] (1) 

where C is a constant combining all the numerical coefficients, p|j(n) are two (s = 1,2) 

polarization tensors and a^{ri), a^{ri), are (Heisenberg) operators for each mode n and for 
each polarization state s. 
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The polarization tensors p|j(n) satisfy the "transverse-traceless" conditions pfjU^ — 0, 
PijS'^^ = and leave independent only two components of hij for each n-mode of the field. 
For a wave travelling in the direction 

n/n — (sin 9 cos </?, sin 9 sin (f, cos 9) the polarization tensors are 

p}j{n) = lilj —niimj, p^j(n) = linij + ljmi, where Ij, rrij are two unit vectors orthogonal to n 
and to each other: Ij = (sin ip, — cos ip, 0), ruj = (cos 9 cos (p , cos 9 sin ip , — sin 9) for 9 < t:/2 
and rrij — —(cos 9 cos </?, cos ^ sin (/?, — sin 9) for ^ > 7r/2. 

The operators a^{ri), a^{r]) are annihilation and creation operators for waves (parti- 
cles) travelling in the direction n. The time evolution of a^{r]), dttiv) is governed by the 
Heisenberg equations of motion for each mode n and for each polarization state s (index 
s is omitted here but will be restored later): dan/drj — —i[an,H], da^/dr] — — i[a+,i7]. 
The Hamiltonian H to be used in these equations has the form H = na^a^ + nal„a_n + 
2(T(77)a+ain + 2(T*(77)a„a_n where the coupling function a{r)) — ia' /2a and ' = d/dr]. The 
solution to the Heisenberg equations of motion can be written as 

aM = ^in(?7)a„(0) + Vn{ni)atM , <{r]) = <(?7)a+(0) + <(?7)a_„(0) (2) 

where an(0), (0)' the initial values of the operators an(r;), a^{rj) taken at some initial 
time long before the coupling became significant and the amplification process has started, 
and the complex functions Um Vn satisfy the equations 

iu'^ — nUn + i{a' / a)v^ , iv'^ — nVn + i{a'/a)u*^ (3) 

where — \vn\'^ = 1 and Mn(0) = 1, Wn(0) = 0. It follows from these equations that 

the function w„ + i>* = //„ obeys the equation ^'^ + (n^ — a" /a)^n — which is precisely 
the equation for classical complex //-amplitude [2] of the gravitational-wave field. Note 
that the solutions m„(?7), Vniji) to Eq. (3) depend only on the absolute value of the vector 

n — {n\ + nl + rig)^/^, not its direction. Also, these solutions are identical for both 
polarizations: they obey the same equations with the same initial conditions. 

The two complex functions «„, f„ restricted by one constraint |MnP — l^nP = 1 can be 
parameterized by the three real functions r„(7^), (pniv)^ ^n{v)'- 
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For each n these functions obey the equations 

r' — {a' /a) cos 2(f), cf)' — —n — {a' /a) sin 2(f) cth2r , e' — —n — {a' /a) sm. 2(f) thr (5) 

which can be used for an exphcit calculation of r„, if a time-dependent scale factor 

a{r]) is given. 

The operators an(0), a,n(0) (Schrodinger operators) satisfy the usual commutation 
relations [an(0), a+ (0)] = (5^(n — m) and the same is true for the evolved operators: 
[a„ (77), 0^1^(77)] = 5^(n — m). By using Eq. (4) the (Bogohubov) transformation (2) can 
be cast in the form 

an(r?) = i?5a„(0)5+i?+ , 0^(7?) = RSa^{<d)S+ (6) 

where 

S{rA) = exp [r (e-2^'^a„(0)a_„(0) - e^"'' al{Q)atM) 

is the unitary two-mode squeeze operator and 

R{e) = exp [-is (a+(0)a„(0) + al„(0)a_„(0))" 

is the unitary rotation operator. The functions r„, 0„, are called squeeze parameter, 
squeeze angle and rotation angle. (For a description of squeezed states see, for example, [4].) 
Equations (2), (6) demonstrate explicitely the inevitable appearance of squeezing in the 
problems of this kind. In this paper we use the presentation based on travelling waves and 
two-mode squeezed states but standing waves and one-mode squeezed states are equally 
good [3]. 

We assume that the quantum state of the field is the vacuum state defined by the 
requirement an(0)|0 >= for each n and for both s. In the Heisenberg picture the state of 
the field docs not change in time but the operators do. The values of a,.^{ri), 0,^(1]) determine 
all the statistical properties of the field at the later times. It follows from Eq. (2) that the 



mean values of an{r]), 0^(77) are zero: < 0|a„(77)|0 >= 0, < 0|a+(77)|0 >= 0, but the mean 
values of the quadratic combinations of a^ir)), a'^{r]) (variances) are not zero: 

< O\an{r])0'm{rj)\^ > = Un{rj)vm{v)^^{T^ + m) 

<0|a+(77)a+ (77)10 > =<(77X(7?)<5^(n + m) (7) 

< 0|an(7?)a+ (77)10 > = Un{v)u*^{r])5^{n - m) 

< 1 a+ (77)0^(77)10 > = vl{ri)vm{v)S^{n - m) 

These relationships (the first two) show explicitely that the waves (modes) with the opposite 
momenta are not independent. On the contrary, they are strongly correlated which is the 
reason for the appearance of standing waves. This fact finds its reflection in the correlation 
functions of the field. 

To simplify the discussion of the correlation functions, we will first ignore the tensorial 
indices in Eq. (1) and consider a scalar field 

/oo 
o?=^n[a„(77)e-^ + a+(77)e--^] . 
-00 

Physically, the field ^1(77, x) may be a scalar variable associated with the density perturbations 
(see Ref. [5] and the third paper in Ref. [3]). The mean value of the field h is zero in every 
spatial point and at every moment of time. The variance of the field is not zero, it can be 
calculated with the help of Eq. (7): 

poo 

< 01/1(77, x)/l(77,x) 10 >= 47r / n^dn(\Un\^ + \Vn\^ + UnVn + u*v*) . 

Jo 

In terms of the squeeze parameters the result can be written as 

roo 

< 0\h{r], x)h{r], x) |0 >= 47r / n^dn{ch2rn + s/i2r„ cos 20„) (8) 

Jo 

(this expression includes the vacuum energy term Att /q°° ri^dn which should be subtracted 
at the end) . The variance of the field does not depend on the spatial coordinate x but does 
depend, in general, on time. The function under the integral in Eq. (8) is usually called the 
power spectrum of the field: P{n) — 'n?{ch2rn + s/i2r„cos20„). The important property of 
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squeezing is that, for a given time, the function P{n) is not a smooth function of n but is 
modulated and contains many zeros or, strictly specking, very deep minima. To see this, one 
can return to Eqs. (5). For late times, that is, well after the completion of the amplification 
process, the function a' /a on the right-hand side of Eqs. (5) can be neglected. (This is 
equivalent to saying that one is considering waves that are well inside the Hubble radius.) 
At these late times, the squeeze parameter is not growing any more and the squeeze angle 
is just (pn — —nr) — 0on- Since r„ ^ 1 for the frequencies of our interest [3], the P{n) can be 
written as P{n) ^ n^e^''" cos^ (nr/ + 0on)- The factor cos^(n?7 + 0on) vanishes for a series of 
values of n; at these frequencies the function P{n) goes to zero. The position of zeros, as a 
function of n, varies with time. The similar conclusions hold for the spatial auto-correlation 
function: 

roo sin Til 

< 0|/i(77,x)/i(77,x-F 1)|0 >= 47r / — — (c/i2r„ s/i2r„ cos 20„)dn . 

«/ lit 

The resulting expression depends on the distance between the points but not on their co- 
ordinates. The power spectrum of this correlation function is also modulated by the same 
factor cos^ {nrj + 0on)- 

We return now to the tensor field (1). There is one combination of the components hij 
which has a special meaning: h{e'^) — hijt^e^ , where 

— (sin 9 cos 0, sin 9 sin 0, cos 9) is an arbitrary unit vector. The h{e^) enters the calculation 
of the CMBR temperature variation seen in the direction (Sachs- Wolfe effect [6]): 

' 2Jo \dr] I 

where w — r)R — r), x'' — e^w, wi — tjr — tje and hij in this formula is a~^{r)) times hij 
introduced in Eq. (1). For a quantized /lij-field, the 6T/T becomes an operator: 
ST 1 /""'i . ^ t 

1 2 Jo J-oo ^^-^ 

+ [<*<+(0)+/?r«i„(0)]e-^'^'=^'-} 

where a^{r}) = (uf^/a)', Pniv) = ("^n/"^)'- The mean value of 6T/T is zero while the variance 
of the expected temperature fluctuations can be written as 
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nT nT 1 /■'"i /""^i 

<Q\ — {e^) — {eH'd>=-C^ dw dw d^ncosinke'^n 
T T 4 Jo Jo J-oo 

X 7r^(n, e*^)/(n, (9) 

where ^ — w — w and 

7r^(n,e'=) ^ (p^,.(n)eVy + (pJ(n)eVy, 
/(n, w, iv) = an{w)al{iv) + + a„(t(;)/3„(iD) + /3*(w)q;*(iD), 

The integration over the variables 6* in Eq. (9) allows one to reduce this formula to 

<0| — (e^) — (e^)|0>=C287r / dw dw n^Wi(n^)f(n,w,w)dn (10) 
T T Jo Jo Jo 

where 

The term 14^1 (n^) depends on the interval between the points but not on the direction 
of sight. Thus, variancies seen in all directions e*^ are the same. They are also position 
independent as for x'^ — e^w + the coordinates Xq of the observer drop out of the final 
result. 

We will now turn to the derivation of the angular correlation function 
< 0|5r/T'(ei)(5T/T(e|)|0 > where e\ and e| are two different unit vectors. The general 
formula for this function can be written as 

ST ST 1 /■'"i /""^'i 

< 0| — (et) — (e2^)|0 > =-C' dw dw d^ncos(n,C) 

11 4: Jo Jo J-oo 

X 7r2(n,et,e^)/(n,«;,w) (11) 

where C,'^ — e\w — e^w and 

7r2(n,e^e^) ^ (pi-(n)elei)(pL(n)e^e-) + (4(n)elei)(pL(n)e^e-). 
The integration over the variables 9 in Eq. (11) reduces this formula to 



01 01 /'^l /""'I 1''-'° 

< 0| — (e^) — (e^)|0 >= C^Stt / dw dw n^W2{nC, coii 5) f{n,w,w)dn 
11 Jo Jo Jo 

(12) 

where C = {w'^ — 'iww cos 5-\-w'^Y^'^, 5 is the angle between the two directions of observation, 

cos 5 = e\e\ + e\e\ + e\e\, and 

« cos 5) = ^ (3 cos^ 5 - 1) (7r/2) V2 (nC) K) 

+ cos5(cos2(5 - l)(nw)(nw)(7r/2)V2(^()-7/2j7/2(nC) 

+ ^(cos^^ - l)\nw)\nw)\n/2Y'\nC)-'''J„2{nC). (13) 
o 

Expression (12) depends only on cos 5 and, hence, the correlation function is rotationally 
symmetric. In the limit cos 5 = 1, the parameter C, goes over into ^ and Eq. (12) coincides 
with Eq. (10). 

Expression (12) gives the angular correlation function in the general and universal form. 
It can be used with arbitrary functions q;„('u;), /5„('u;), that is, it is applicable for arbitrary 
(not necessarily inflationary) cosmological models generating squeezed gravitational waves. 
The remaining integrations in Eq. (12) assign concrete numerical values to the correlations 
attributed to different separation angles 5, but they do not change the general angular pat- 
tern represented by the function W2{nC,, cos 5). Consistency with the data of the COBE-type 
observations may lead to the determination of the functions an{w) , (3n{w) and, eventually, 
to the knowledge of the expansion rate of the early universe. This will be a subject of a 
separate discussion. The implications of the COBE observations for inflationary models are 
under active analysis (see, for instance, a recent paper [7] and references therein) . Some new 
results based on the correlation function (12), (13) have been derived in [8]. 
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